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Предлагается подход к определению правил функционирования клеточных автоматов, отличающийся совместным использова- 
нием четырех пространственных характеристик и позволяющий увеличит адекватность моделирования. Приведен сопостави- 
тельный анализ предложенного подхода для решения задачи анализа динамики и прогнозирования ландшафтного покрова на 
примере природно-территориального комплекса Португалии. 


Введение 

В связи с постоянно растущими техническими 
возможностями вычислительной техники, аппа- 
ратного, алгоритмического и программного обес- 
печения космического и наземного сегментов со- 
временных систем дистанционного зондирования 
Земли (ДЗЗ) увеличиваются возможности решения 
сложных задач анализа и прогнозирования разви- 
тия явлений различных природно-территориаль- 
ных комплексов (ПТК). Интересным и сложным 
для изучения динамики и прогнозирования пред- 
ставляется такой компонент ПТК как земная по- 
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верхность, и, в частности, ландшафтный покров, 
отличающийся существенной пространственной 
распределенностью и динамичностью изменения 
своих границ [1]. 

Наиболее часто анализ динамики и прогнозиро- 
вание изменений ландшафтного покрова выполня- 
ется с использованием детерминированных и сто- 
хастических моделей [2]. Однако на сегодня ш ний 
день существует лишь несколько удачных примеров 
реализации таких моделей, использование которых 
оправдано при практическом решении задач изуче- 
ния динамики и прогнозирования изменений ланд- 
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шафтного покрова [3-7]. Кроме того, существует 
проблема выбора модели, которая наиболее подхо- 
дит для исследования и позволит обеспечить аде- 
кватное прогнозирование изменений ландшафтно- 
го покрова в каждом конкретном случае. 

В данной работе проводится сравнительный 
анализ современных средств для решения задач 
прогнозирования изменений ландшафтного по- 
крова, включая идентификацию проблем, препят- 
ствующих широкому практическому использова- 
нию таких средств. Предлагается подход, отличаю- 
щийся совместным использованием четырех про- 
странственных характеристик и позволяющий уве- 
личить адекватность моделирования. Обсуждаются 
результаты применения предложенного подхода и 
других средств для решения задачи анализа дина- 
мики и долговременного прогнозирования ланд- 
шафтного покрова ПТК Синтра-Кашкаиш (Порту- 
галия) и даются рекомендации по выбору моделей 
изменения ландшафтного покрова, а также путей 
их использования при решении практических за- 
дач. 

Существующие подходы к моделированию 

изменений ландшафтного покрова 

Математический аппарат моделей, используе- 
мых в задачах анализа динамики и прогнозирова- 
ния изменений ландшафтного покрова, в дискрет- 
ном пространстве состояний представляют с ис- 
пользованием разностных уравнений в матричной 
форме записи [2]: 

и (+1 =Рхи„ (1) 

где и ( =[и 1 ,я 2 ,...,я т ] т - вектор-столбец абсолютных 
вероятностей того, что система перейдет в одно из 
т состояний к моменту времени (/+1), Т - знак 
транспонирования, Р - квадратная матрица поряд- 
ка т, каждый элемент рд которой является вероят- 
ностью перехода элемента из состояния 0, в со- 
стояние 0, за временной интервал, начиная с мо- 
мента времени /до (/+1). 

Принимая во внимание сложность анализа эко- 
систем в целом, и ландшафтного покрова, в част- 
ности, отмечают несовершенство имеющегося ма- 
тематического аппарата, алгоритмического и про- 
граммного обеспечения существующих моделей, 
не позволяющих в достаточной степени эффектив- 
но решать задачи анализа и прогнозирования изме- 
нений ландшафтного покрова. Подавляющее боль- 
шинство моделей, реализованных в качестве моду- 
лей и подсистем в различных программных систе- 
мах, представляют собой исследовательские вер- 
сии, не подвергавшиеся тщательному анализу. По- 
этому внимание исследователей в области анализа 
динамики ландшафтного покрова привлекают нес- 
колько наиболее законченных программных реа- 
лизаций моделей: СШЕ ( Сопѵетоп о/ Ьапсі ІІж апё 
Из Е//есІз), Ь ІІСА5 ( ЬапсІ ІІзе Скап§е Апаіуш ЕузГет), 
СЕМ (Сепегаі Есоауаіет МосІеІ), РЬМ (РаІихепГ ЬапсІ - 
жаре МосІеІ) и НОМ ( ІІгЪап 6го\ѵік МосІеІ) [8]. 


Среди этих реализаций моделей следует выде- 
лить ГГСМ, предназначенную для исследования 
роста урбанизированных территорий и предста- 
вленную в программном модуле САМагкоѵ в мно- 
гофункциональной растровой геоинформацион- 
ной системе (ГИС) Мгш32 (Сіагк ЬаЬз, Сіагк Ші- 
ѵегзііу, США) [9]. Эта модель (будем называть ее по 
имени реализующего ее программного модуля - 
СА Магкоѵ) основана на использовании клеточных 
автоматов (КА), позволяющих учитывать про- 
странственные особенности ландшафтно-классо- 
вой структуры при моделировании [10]. 

Другая модель, представленная в ГИС Ыгі$і32 
Кііітапіаго и реализованная в модуле Сеотосі (назо- 
вем эту модель по аналогии - Сеотосі), в некоторой 
степени представляет собой упрощенную версию 
модели СА_Магкоѵ. Модель Сеотосі позволяет про- 
водить при моделировании замену только одного 
ландшафтного типа на другой и не использует в 
своей основе КА, а пространственные особенности 
ландшафтных типов учитываются через опцио- 
нальное применение фильтра в скользящем окне. 

Сложность моделируемых процессов, отсут- 
ствие достаточной априорной информации, ис- 
пользование часто слишком простых и единых на 
всем моделируемом изображении правил функцио- 
нирования КА (или параметров фильтрации) не по- 
зволяют учитывать многие пространственные осо- 
бенности исследуемого ландшафта, что в конечном 
итоге может вести к снижению адекватности моде- 
лирования и сложностям практического примене- 
ния таких моделей как СА Магкоѵ и Сеотосі. 

Подход к определению правил функционирования 

клеточных автоматов 

Одной из наиболее сложных задач при проведе- 
нии моделирования с помощью КА является опре- 
деление правил их функционирования. Так, в [7] 
правила функционирования КА предлагается 
определять в каждой точке исследуемой террито- 
рии на основе результирующей вероятности заме- 
щения /-го на у'-ый тип ландшафта: 

РГ =/(РіГ’ РіГ> РГ)> ( 2 ) 

где рр" - вероятностная характеристика, определя- 
емая на основе значений стохастической матрицы 
Р, />/ оп - вероятность, рассчитываемая на основе 
априорной информации, а -вероятность, значе- 
ние которой базируется на одной пространствен- 
ной характеристике — критерии насыщенности 
[11]. Применение нескольких пространственных 
характеристик для определения р- п в (2), позво- 
ляющих более полно учесть особенности ландшаф- 
та, как правило, затруднено вследствие различной 
информативности этих характеристик и суще- 
ственно различающихся диапазоном их значений. 
В этой связи предлагается подход к вероятностно- 
му определению правил функционирования КА, 
основанный на использовании четырех наиболее 
информативных пространственных характеристи- 
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ках из [12], совместное применение которых при 
определении Р;"" обеспечивается единым диапазо- 
ном изменения их значений — [0; 1] , а также за счет 
эквивалентных весовых коэффициентов. 

В данном случае расчет пространственных осо- 
бенностей ландшафта производится в скользящем 
по моделируемому изображению окне порядка й с 
количеством элементов в нем (2сІ+1)х(2с1+\) (при 
сІ=\ размер окна 3x3). Рассмотрим более подробно 
расчет четырех характеристик с учетом того, что 
каждый компонент этих характеристик рассчиты- 
вается в выбранном скользящем окне порядка й. 


Первая характеристика, показывающая частоту 
встречаемости к-то ландшафтного типа, рассчиты- 
вается следующим образом: 


Ъ п ш 
’ = 2=1 
ЧВ(*,І) 


( 3 ) 


где п ки - число элементов к-то типа ландшафта у-го 
участка в рассматриваемой окрестности /- ой точки 
изображения, п, - общее число ячеек в рассматри- 
ваемой окрестности, /и, - общее число всех участков 
в окрестности. Под участком будем понимать набор 
элементов одного типа ландшафта, окруженных эл- 
ементами типов ландшафта, отличными от рассма- 
триваемого. Если исследуемая окрестность будет со- 
стоять из одного типа ландшафта, то С чвад = I . 

Вторая характеристика, показывающая степень 
фрагментации А;-го типа ландшафта, рассчитывает- 
ся следующим образом: 


где т кі — число участков к-то типа ландшафта в 
окрестности /-ой точки. Если исследуемая окрест- 
ность будет содержать два участка для разных ти- 
пов ландшафта, то С ФЛ , ■,=(), 5. 

Третья характеристика, показывающая среднее 
фрактальное измерение периметров к-то типа 
ландшафта, рассчитывается как: 


С, 


X 

7=1 


21п(0,25х^, . .) 

Нп к .,) 


СФЩк,і) 


— 1, 


Ш. 


(5) 


где р кц - периметр у-го участка А;-го типа ландшаф- 
та в окрестности /-ой точки. Если все участки к-то 
типа ландшафта являются простейшими (напри- 
мер, квадрат), то С сфи( *, 0 =0. 

Четвертая характеристика, показывающая 
среднее расстояние между участками к-то ланд- 
шафтного типа, рассчитывается как: 

ш і 

2Хм 

С„„„ = (6) 

где Іг ки - число элементов от у-го участка к-то типа 
ландшафта до ближайшего соседнего элемента 


участка к-то типа в окрестности /-ой точки, сі - по- 
рядок скользящего окна. Если участки одного 
класса отдалены друг от друга на максимальное 
расстояние, то С СБР(М =1. 

Для расчета вероятностей замещения в тот или 
иной тип ландшафта для каждой /-ой точки изо- 
бражения в соответствии с выражениями (3)-(6) 
предлагается использовать два вида характеристик. 
Первый вид рассчитывается для текущей анализи- 
руемой окрестности (локальные характеристики), 
второй - для всего изображения в целом (глобаль- 
ные характеристики). 


Локальные характеристики для каждого из ти- 
пов ландшафтного покрова рассчитываются в теку- 
щей анализируемой окрестности с использованием 
выражений (3)-(6), и могут быть представлены в 
виде вектора 


С лок(к) 


с 

'“ / ЧВ (к,і) 

с 

ФОМ) 

г 

^ФЩк,і) 

с 

^СБРОМ) 


(7) 


где к - тип ландшафтного покрова, / - текущая 
анализируемая точка изображения. 

Глобальные характеристики могут быть пред- 
ставлены в виде вектора 


Сгд(к) 


г 

'-'ЧВ(А-) 

с 

^Ф(к) 

с 

^ФИЦ) 

с 

'-'СЕРЦ) 


(8) 


где каждый его компонент рассчитывается для к-то ти- 
па ландшафтного покрова путем перемещения сколь- 
зящего окна по всему изображению и подсчета средне- 
го значения для каждой из характеристик (3) — (6). 

Сравнивая значения векторов локальных и гло- 
бальных характеристик для различных типов ланд- 
шафта можно определить, насколько значения 
пространственных характеристик в скользящем 
окне схожи со значением пространственных харак- 
теристик на исследуемой территории в целом. Это 
позволит выявить вероятность принадлежности 
центрального элемента анализируемой окрестно- 
сти к к - му типу ландшафта. Учитывая это, про- 
странственную вероятностную составляющую 
предлагается рассчитывать как: 

Р цр — |1— й?(Р С док(і)> Рсгл№))1’ (9) 

где сІ(Р С:ж щ, Р Сгл№) ) - Евклидово расстояние между 
векторами (7) и (8). 

Использование в выражении (2) вероятности 
из (9) позволит более точно учитывать простран- 
ственные особенности ландшафта и адекватно рас- 
считывать правила функционирования КА в каж- 
дой точке исследуемой территории. 
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Прогнозирование изменений ландшафтного 
покрова с использованием различных моделей 

Территория исследования и используемые исход- 
ные данные. Решение задачи анализа динамики и 
прогнозирования выполнялось для ПТК Синтра- 
Кашкаиш площадью 416 км 2 (рис. 1), расположен- 
ного на прибрежной территории к северо-западу от 
г. Лиссабона (Португалия). Эта территория харак- 
теризуется высоким уровнем плотности населения 
и его существенной динамики прироста (в период с 
1990 по 2000 гг. рост составил около 30 %), что су- 
щественно повлияло на изменение границ и пло- 
щадей ландшафтных классов, на места и условия 
проживания большого количества людей. 



Рис. 1 . Изучаемая территория - Синтра-Кашкаиш (Португа- 
лия) 


Для решения задачи анализа динамики и прог- 
нозирования изменений ландшафтного покрова с 
помощью моделирования обычно используют нес- 
колько ландшафтных карт, полученных на основе 
космических снимков с систем ДЗЗ или путем кар- 
тирования на местности. При этом для расчета па- 
раметров моделирования используют карты по- 
верхности для моментов времени (7—1) и 7, а соб- 
ственно моделирование выполняют на основе кар- 
ты для момента времени 7, принимая ее в качестве 
исходной. Кроме того, используют карту для мо- 
мента времени (7+1) для оценки результатов моде- 
лирования и адекватности созданной модели. 

В качестве исходных данных для моделирования 
изменений ландшафтного покрова использовались 
разновременные космические снимки системы ДЗЗ 
Ьапсізаі ТМ (получен 08.02.1994 г.) и ЬапсЕаІ ЕТМ+ 
(получены 14.03.1989 г. и 08.04.2001 г.), что позволя- 
ет осуществить моделирование с использованием 
карт, полученных по космическим снимкам 1989 и 
1994 гг., а оценку адекватности моделирования вы- 
полнить с использованием карты, полученной по 
космическому снимку 2001 г. 

Для каждого из космических снимков, для та- 
ких ландшафтных классов как заселенные и неза- 
селенные территории, по имеющимся априорным 
данным наземных наблюдений были построены 
обучающие выборки. Снимки были сегментирова- 
ны с использованием алгоритма иерархической 
кластеризации, представленного в программном 
обеспечении системы еСо§пШоп 3.0. Для этого ис- 


пользовались все каналы (диапазоны съемки) кос- 
мических снимков ДЗЗ Ьапскаі (за исключением 
теплового), а также применена пороговая характе- 
ристика, матрица смежности и текстурная характе- 
ристика «однородность», что позволило с доста- 
точно высокой точностью относить каждый из сег- 
ментов либо к классу «заселенная территория», ли- 
бо - «не заселенная». 

При моделировании изменений ландшафтного 
покрова обычно используют имеющуюся априор- 
ную информацию о влиянии тех или иных факто- 
ров на развитие (угнетение) ландшафтных типов на 
исследуемой территории. Формализация этой ин- 
формации осуществляется с помощью функций 
пространственного анализа современных ГИС пу- 
тем построения вероятностных карт соответствия 

[13] . В данной работе для построения таких карт на 
исследуемую территорию были использованы ци- 
фровые карты расстояний до автодорог и населен- 
ных центров, а также другая дополнительная ка- 
дастровая информация, полученная из Португаль- 
ского географического института, а также цифро- 
вые карты уклонов, полученные со спутника НАСА 
(ЗНиІіІе Каёаг Торо§гарНу МШіоп). 

Критерии оценки точности. На сегодняшний 
день серьезным препятствием при анализе дина- 
мики и прогнозировании изменений ландшафтно- 
го покрова является до конца нерешенная пробле- 
ма оценки адекватности используемой модели и 
проблема оценки точности полученных результа- 
тов моделирования (прогноза). Эти проблемы свя- 
заны с тем, что для оценки точности полученных 
результатов моделирования необходимо иметь ре- 
альные (не модельные) данные на соответствую- 
щий момент времени. Однако даже при наличии 
таких данных полученная количественная оценка 
точности не всегда может адекватно отражать суть 
процессов, происходящих на земной поверхности 

[14] . Следует особо отметить визуальный и числен- 
ный способы оценки результатов моделирования, 
которые в некоторых случаях могут давать различ- 
ные результаты. 

Для численных способов оценки точности, как 
при классификации аэрокосмических изображе- 
ний, так и при моделировании изменений процес- 
сов ландшафтного покрова применяют самые раз- 
личные критерии. Среди наиболее используемых 
исследователями следует отметить такие критерии 
как каппа индекс согласия (КИС), оѵегаіі ассигасу, 
Кпо, Кіосаііоп и КНШо [15]. Каждый из таких крите- 
риев обеспечивает оценивание точности результа- 
тов по изображениям или тематическим картам с 
учетом своих особенностей. Например, Кіосаііоп 
оценивает особенности местоположения элемен- 
тов ландшафтного покрова сравниваемых темати- 
ческих карт, КНШо обеспечивает количественное 
сравнение элементов ландшафтного покрова этих 
карт. При этом КИС рассчитывается как произве- 
дение Кіосаііоп и КНШо, что позволяет учитывать и 
пространственные и количественные особенности 
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Рис. 2. 
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Фрагменты исходных тематических карт: 
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а) 1989 г., 6) 1994 г., в) 2001г. 


сравниваемых фрагментов карт. В этой связи 
облегчить задачу оценки точности результатов мо- 
делирования и сделать ее более объективной следу- 
ет с помощью совместного использования указан- 
ных выше критериев точности. 

Построение вероятностной карты соответствия. 
Как отмечалось выше, для моделирования целесо- 
образно использовать имеющуюся априорную ин- 
формацию, которая может быть учтена с помощью 
вероятностных карт соответствия. Построение та- 
ких карт в этом случае следует выполнить лишь для 
двух типов ландшафтного покрова, для чего были 
использованы цифровые карты расстояний до ав- 
тодорог и населенных центров, а также цифровые 
карты уклонов. Кроме того, на основе растровых 
карт на 1989 г. и 1994 г. (рис. 2, а, б) были построены 
дополнительные вероятностные карты соответ- 
ствия допуская то, что территория, заселенная в 
1989 и 1994 гг. будет заселена и в будущем. 



Рис. 3. Общий вид интегрированной вероятностной карты 
соответствия 

В заключение из вероятностных карт соответ- 
ствия путем перемножения соответствующих веро- 
ятностных значений и нормирования была сфор- 
мирована интегрированная вероятностная карта 
соответствия, содержащая в каждой точке исследу- 
емой территории информацию о вероятности ее за- 
селения (рис. 3). 


Используемые модели. Выше отмечалась перс- 
пективность практического применения моделей 
СА_Мсігкоѵ и ОеошосІ, реализованных в ГИС Ыті 
Кііітапіаго. Учитывая ряд отмеченных недостатков, 
присущих моделям этого типа, была предложена и 
реализована модель, условно названная СА АсІѵап- 
сесі. Вычислительная часть этой модели реализова- 
на в среде Місгозой Ѵізиаі С++ 6.0, компилятор ко- 
торой обеспечивает высокую степень оптимизации 
двоичного кода и гибкую работу с памятью, что по- 
зволило обеспечить высокую вычислительную эф- 
фективность. Кроме того, в СА_АсІѵапсеб реализо- 
ван подход к определению правил функционирова- 
ния клеточных автоматов, предложенный в данной 
работе и описанный выше. Все это позволило при 
проведении исследований, направленных на изу- 
чение динамики и на прогнозирование изменений 
ландшафтного покрова территории Синтра-Каш- 
каиш, использовать три упомянутых модели - 
СА_Магкоѵ, ОеотоА и СА_Абѵапсеё и оценить пред- 
ложенный подход на практике. 

Необходимым этапом, предваряющим практи- 
ческие исследования с использованием моделей, яв- 
ляется этап их калибрования. На этом этапе опреде- 
ляются по возможности оптимальные параметры 
модели - размер скользящего окна, количество ите- 
раций модели, соответствующее необходимому ре- 
альному временному промежутку и т. п. В данном 
случае калибрование моделей с выбором наилучших 
параметров осуществлялось следующим образом. 
Параметры моделей рассчитывались по имеющимся 
растровым картам для 1989 г. (момент времени (1-1)) 
и для 1994 г. (момент времени 1). Затем в качестве ба- 
зового (исходного) изображения использовалась 
карта на момент времени (1-1), а не на момент вре- 
мени 1, как обычно при прогнозировании. Это по- 
зволило проводить сравнение полученных прогноз- 
ных карт на момент времени 1с реальными данными 
на этот же момент времени и определить наилучшие 
параметры моделей в каждом случае. 

Результаты исследований. Учитывая целесооб- 
разность оценивания точности моделей с исполь- 
зованием различных критериев, при проведении 


84 




Технические науки 


исследований в данной работе совместно примене- 
ны несколько статистических показателей (см. та- 
блицу). Численная оценка величины изменений, 
произошедших за период с 1994 до 2001 гг. прово- 
дилась путем сравнения соответствующих темати- 
ческих карт на эти годы. 


Таблица. Результаты оценки точности прогнозных карт на 
2001 г. для разных моделей, % 


Критерий точности 

Модели | 

СА_Магкоѵ 
окно 5x5 

СА_Асіѵапсесі 
окно 3x3 

СеотосІ 
окно 5x5 

« Ноль -модель » 
1994 г. 

ОѵегаІІ ассигасу 

83,59 

83,95 

82,96 

80,00 

КИС 

64,82 

66,28 

64,02 

57,35 

Кпо 

67,17 

67,90 

65,92 

59,99 

Кіосаііоп 

75,27 

68,57 

68,31 

64,57 

КЬЩо 

86,12 

96,66 

93,72 

88,82 


Результаты в таблице показывают, что более вы- 
сокая адекватность моделей изменения ландшаф- 
тного покрова достигается в случаях, если происхо- 
дящие (моделируемые) изменения достаточно су- 
щественны и носят тенденциозный характер. Ре- 
зультаты, представленные в столбце 5 таблицы по- 
казывают относительно высокие показатели для 
«ноль-модели» (зарубежные исследователи исполь- 
зуют понятие «Шй тоёеі») означающие, что грани- 
цы и площадь типа «заселенные территории» с 1994 
по 2001 гг. изменились не очень существенно, что, 
в свою очередь, позволяет отметить сложность 
оценки любых моделей с использованием имею- 
щихся в данном случае разновременных данных. 
Несмотря на отмеченную сложность, практически 
все полученные оценки точности для модели 
СА_АсІѵапсесІ (выделено жирным шрифтом) нес- 
колько выше оценок для двух других моделей, что 
позволяет сделать вывод о целесообразности и до- 
статочно высокой эффективности применения 
подхода к моделированию с помощью КА с ис- 
пользованием пространственных характеристик 
(3) — (6). Кроме этого, можно отметить перспектив- 
ность проведения более масштабных дополнитель- 


ных исследований, направленных на анализ ин- 
формативности пространственных характеристик 
при определении правил функционирования КА. 

Долговременный прогноз. Задачи анализа дина- 
мики ландшафтного покрова предполагают прове- 
дение как кратковременного (2-5 лет), так и долго- 
временного (5 и более лет) прогнозов. Проведение 
долговременного прогноза позволяет оценить ра- 
ботоспособность и эффективность оцениваемых 
моделей, т. к. при практическом применении моде- 
ли, пригодные для краткосрочного прогноза, не 
пригодны для долгосрочного, и наоборот. Несмо- 
тря на существенную непредсказуемость развития 
ландшафтного покрова в длительном периоде, дол- 
говременный прогноз с использованием работос- 
пособной модели может быть полезен и позволит 
эксперту выявить и оценить существующие тен- 
денции ландшафтного покрова и позволит ему 
предпринять действия, направленные на снижение 
вероятности возникновения кризисных ситуаций 
на исследуемой территории в будущем. 

В данной работе на моделях СА_Магкоѵ, СА АсІ- 
ѵапсесі и Сеотод был проведен ряд экспериментов 
по проведению долговременного прогнозирования 
до 2025 г, некоторые результаты которых предста- 
влены на рис. 4. Для возможности сравнительного 
анализа результатов моделей использовались иден- 
тичные исходные данные. Имеющиеся растровые 
карты заселенных территорий и кадастровая ин- 
формация позволила представить исследуемую 
территорию в виде 26 административно-террито- 
риальных единиц в каждой из которых была оцене- 
на значимость ожидаемых изменений. 

В связи с тем, что прирост заселенных террито- 
рий (рис. 4, а) не носит поступательного характера, 
полученные результаты позволяют сделать вывод о 
недостаточной применимости модели СА_Магкоѵ 
для выполнения долгосрочного прогноза. В то же 
время использование СА_АсІѵстсесІ и Оеотосі дают в 
данном случае более реалистичные результаты (за 
исключением одной из территорий, где модель Се- 


Динамика 
заселения (%) 


[]< 15 

(И (-15- 'ІО! 

■ (-Ю-5) 

1 I (~5 ~ 5] 

I | (5 - Ю] 

Л (10- 15] 
И (15 -20] 
Щ (20 - 25] 
^■>50 



Рис. 4. Карты прогноза роста заселенных территорий региона Синтра-Кашкаиш (Португалия) к 2025 г., полученные с исполь- 
зованием моделей: а) СА_Магкоѵ, б) СА_Ас/ѵапсес/, в) беотосі 
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отосі «прогнозирует» неоправданно высокий рост 
плотности заселения). Учитывая то, что Оеотос! 
обеспечивает моделирование поведения только од- 
ного ландшафтного класса, при сложной ландшаф- 
тно-классовой структуре исходных данных исполь- 
зование этой модели может быть не оправдано. 

Заключение 

Задачи анализа динамики ландшафтного по- 
крова, краткосрочного и долгосрочного прогнози- 
рования роста заселенных территорий даже сегод- 
ня являются крайне сложными, а существующие 
алгоритмические и программные средства не всег- 
да позволяют исследователям решать такие задачи 
с достаточной эффективностью. Данная работа бы- 
ла направлена на поиск более эффективных спосо- 
бов и средств решения задач анализа динамики и 
прогнозирования изменений ландшафтного по- 
крова, в ходе которой были получены следующие 
основные результаты: 
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1. Предложен и программно реализован подход к 
определению правил функционирования КА, 
который при моделировании позволяет более 
точно учитывать пространственные особенно- 
сти ландшафта, чем в моделях СА_Магкоѵ и Се- 
отосі, увеличивая адекватность моделирования. 

2. Проведен сравнительный анализ моделей 
СА_Магкоѵ, СА_ААѵапсесІ и Оеотос і при решении 
практической задачи анализа динамики, крат- 
косрочного и долгосрочного прогнозирования 
роста заселенных территорий региона Синтра- 
Кашкаиш (Португалия), который показал перс- 
пективность предложенного подхода, особенно 
при долгосрочном прогнозировании в условиях 
сложной ландшафтно-классовой структуры. 

3. Сделан вывод о том, что для решения задачи 
долговременного прогнозирования в условиях 
сложной ландшафтно-классовой структуры не- 
целесообразно использование моделей 
СА Магкоѵ и Сеотоё. 
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